
function [dfdt] = waves_benchmark(t,y)
global L;

%
dfdt = -L*y;
